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Abstract 

Evolution of extended data is considered in various flow problems, using Nambu brackets 
as a tool applicable to all cases. Extra dimensions, N-brackets, and extended structures 
are first employed to linearize the Euler equations. N-bracket induced evolutions of strings, 
membranes, and other d-branes are then discussed in detail. 



ANL-HEP-CP-03-058 A talk given at the International Conference on Nonlinear Evolution Equations and 
Applications, Northwestern University, 12 June 2003, based on work with David Fairlie, University of Durham, and 
Cosmas Zachos, Argonne National Laboratory. Thanks to Ernest Hemingway for the mojito-inspired title. 
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This talk consists of two more or less independent parts. The first part is a discussion of simple 
Euler fluid flow, in somewhat unconventional terms. The second part is a consideration of the 
motion of continuous collections of points (i.e. branes) in response to a given ambient flow field. 
Although treated independently, both considerations involve some common mathematical language: 
the use of Nambu brackets and extended data structures. 



1 Linearizing the Euler-Monge equations 



General results in n dimensions Encode the local velocity field for n-dimensional flow u (x, t) 
into a non-local structure, after dimension-doubling. Define 

gafcUfc(q,t) _ 

dqi--- dqn 6 (g^ - Xk) ( I x 



X det 



Uk (x, t, a) = 



^ {e {qn - Xn) e-"^"'^^''^) J 



:i) 



exclude kth row 
and kth column 



where e [s) = ±| for s ^ . The local velocity data is given by the extra-dimension boundary 
values 

Uk (x,i(:) = lim Uk (x,t,a) . 



Then with 

we have 

if and only if 



Mn [U] 



d_ 

dt 



d 



^ dxi 



Hn = 



d_ 

dt 



dxjdtti 



Mn[u\ Ui (X,t) = 



7^„f/fc(x,^,a) = . 

The latter equations are solved by the heat kernel method, as given formally by 

Uk (x, t, a) = e ^ Uk (x, t = 0, a) . 



(2) 
(3) 

(4) 
(5) 

(6) 



Expansion and evaluation of the RHS yields the formal time power series solution in terms of the 
initial data. Restricting this series to the extra dimension boundary, a = 0, we obtain the time 
power series solution for Uk (x, t) in terms of Uk (x, t = 0) and its derivatives. 

Dimension-doubling here is formally similar to point-splitting in Schrodinger quantum mechanics 
pp. But it is nonetheless a bona fide doubling of dimensions, in the same sense as the Fourier 
transform of point-split wave-function bilinears in quantum mechanics yields a Wigner function on 
a phase-space whose dimension is twice that of the original wave-function configuration space. 
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Nambu brackets Classically, N-brackets are defined as Jacobians. 

These multi-linear generalizations of the Poisson bracket arise naturally in the previous construction, 
since 



S (qk - Xk) X det (■ ■ ■ )exciude kth row = {s {qi - Xi) 6^^"% ■ ■ ■ , E {qk - Xk) , ■ ■ ■ ,e{qn-Xn) e""""}„^...„ , 

and kth column 

(8) 

with the exponential missing from the kth entry in the bracket. 

In the extra-dimension-boundary limit, the N-bracket becomes a Dirac-delta equating all com- 
ponents of X to those of q. 

hm {s (gi - x^) e'^^^\ . . . , ^ (g, _ a;,) , ■ ■ ■ , e (g„ - a;„) e""""}^^...^^ = ^("^ (q - x) . (9) 

This ensures that U collapses down to u on the a ^ boundary. 

Note the Nambu bracket is a derivation, and is linear in each argument. So another way to say 
the above is 

S {qk - Xk) 6"'="'= det (■ ■ ■ )exciude fcth row = {e {qi - Xi) e"^"S ■■■ ,e{qn- Xn) e""""}„^...„ 

and fcth column 

-ake{qk-Xk) e"'="'= {e (gi - Xi) e''i"\ ■ ■ ■ , m^, ■ ■ ■ , £ (g„ - x„) e«"""}^^...^^ . (10) 
A trivial example in one-dimension As the simplest, local example, consider 

^au(x,t) _ ]^ 

U{x,t,a) = , u{x,t) = —\ia{l + all {x,t,a)) = limU {x,t,a) . (11) 

a a o,~*o 



Then 



if and only if 



Moreover, 



g-U{x,t,a) = -^U{x,t,a) (12) 
d d 

—u{x,t)=u{x,t)—u{x,t). (13) 



U{x,t,a) = e^9^ U{x,t = 0,a) (14) 
immediately gives the time power series solution from initial data 

2 / pau{x) _ 1 \ °o p (ji 

[—^) . (15) 

So, for example, an initial linear profile, u{x) = a + Px, evolves into u (x, t) = (a + Px) / (1 — Pt). 
For a given interval in u, this steepens as time increases, for /3 > 0, and becomes vertical at tbreak = 
1//3. Other simple breaking wave examples can be found in [T] (or at http:/ /curtright. com/waves. html) ) . 



Two-dimensions and some non-locality Define the extended data structures {Ui = U, U2 



V, etc] 



U (x, y, t, a,b)= [dreiy- r) e''<-^''^'^+^-(-^r,t) du (x,r,t) 



dr 



V{x,y,t,a,b) = [dqe{x-q) e-"(^'^'*)+H'?.?/,t) ^^^^'^'^^ . 

J oq 



(16) 



Then 

if and only if 



= = (17) 



M2U = = M2V . (18) 
Again, note the extra-dimension boundary limits 

u{x,y,t)= Urn U {x,y,t,a,b) , v{x,y,t)= \im V {x,y,t, a, b) . (19) 

The demonstration of the equivalence of (jl7|) and (jl8|) goes as follows. Direct calculation gives 

n2U (x, y, t, a, b) = e'^"(-'S/'*)+^-(-'S/.*)>l2M (x, y, t) 

+ b I dre{y-r) e-(-''-'*)+^-(-.-.*) ( ^!^i^2ll^M2V {x, r, t) - ^^i^M^^^^ (3;, r, t) 
J \ or or 

Similarly for V. So then Ai2U = = Ai2V clearly implies 7^2^^ = 0. The converse follows by 
using the obvious limit 

\im T-C2U {x,y,t,a,b) = M.2U {x,y,t) , (20) 

a,b—>0 

along with the similar relation for V. qed 

Three-dimensions and Poisson brackets (Ui = U, f/2 = V, f/3 = W, etc.) 

HsU = HsV = HsW = (21) 

if and only if 

M3U = M3V = Msw = , (22) 

where 



U (x, y, z, t, a, b,c) = / dr e{y — r) e 



au+hv+cw l-^' ^' ^' ^) 



-// 



dr 

drds e{y~r) e{z-s) e""+'»'+^'^ {«, w}^^ (x, r, s, t) , (23) 



V (x, y, z, t, a, b,c)= I dse{z~ s) e'^"+'»'+^'" i^^ V^s^t) 

OS 

au+bv+cw 



a 1 1 dqds e{x-q)e{z-s) 6-+"^+^- {v, u}^^ (g, y, s, t) , (24) 



W {x,y, z,t,a,b,c) = J dqe{x — q)e' 



au+bv+cw 

dw {q,y,z,t) 
dq 



dqdr e{x-q) e{y-r) e'^"+^''+^^ {w, v}^^ {q, r, z, t) , (25) 



I qr 

and where the Poisson bracket is given as usual by 

dudv dudv 

Once again the proof of the equivalence of ^2\\ and ()22j) is given by direct calculation. Higher 
dimensions lead to higher rank Nambu brackets, as expressed by the general n result quoted above. 



2 Nambu evolution of d-branes 

Consider flow-driven motion in a local velocity field u [x] . For points moving along with the flow^ 

dx' = u' [x] dr . (27) 

This deceptively simple description is general enough to incorporate the phase-space evolution of 
any Hamiltonian dynamical system, upon identifying pairs of variables as canonically conjugate, 
say {x^,x^^^) — > {x^p'). 

Is there an elegant way to describe the motion of an extended, continuous collection of points 
moving with the flow? For example, consider (vortex) strings, (surface interface) membranes, or 
higher dimensional manifolds of points immersed in fluids flowing in higher dimensions. The latter 
higher dimensional generalization is most important for dynamics in phase-space, which for a system 
with degrees of freedom is 2N dimensional. That is to say, if we parameterize the extended set 
of points for a given r as 

X {T,ai,a2, - ■ ■ ,Oid) , (28) 

then we have a d-dimensional data-brane, or "d-brane" for short. What is the most economical 
way to describe the collective evolution of this extended set of data? 

A suggested answer to both these rhetorical questions is as follows. If the ambient velocity field 
is given by a Nambu bracket, as 

u'[x] = {x\h,--- Jn-l} , (29) 

where the Is are "flow invariants," then we have satisfied a sufficient condition to obtain an exact 
n-form, Q = dA, that provides a minimal coupling of the d-brane to the velocity field, for the case 
d = n — 2. The individual points are governed by Nambu dynamics [S] in this case, with 

dx^ 

— = {xMi,---,4_i} . (30) 

^There is one caveat to keep in mind here. It is tempting to think of r as conventional time t, defined universaUy. 
But this need not be so. In fact, in some situations it can be misleading to think of r this way. For streamline fluid 
flow, e.g., T could be allowed to vary from one streamline to the next such that evolution along each streamline is 
governed by that streamline's own individual clock. In this sense, t is more like proper time in relativity. In the 
context of Nambu mechanics in phase-space, r is generally related to t by a dynamical scale factor, and therefore 
cannot be simply identified with universal time across all dynamical sectors [2]. That is to say, t may be trajectory 
dependent. For this and other reasons to become apparent, we will call t "Nambu time." 
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Examples of such dynamical systems in phase-space are abundantly provided by superintegrable 
systems. An N degree-of-freedom system, with phase-space of dimension 2N, is maximally super- 
integrable and described as above when n = 2N with a corresponding set of 2N — 1 invariants, /j. 
In this 2N phase-space, there can be no more than 2N — 1 invariants. These ideas are extensively 
discussed in 0. 

As already noted, the preferred continuous collection of points, i.e. the minimally coupled d- 
brane, is of dimension d = n — 2 when the velocity field is given by the above Nambu bracket. The 
explicit n-form that provides the minimal coupling of the brane to the velocity field is then 

Vt = ei^...i„ {dx^^ - vj-^ [x] dr) A ■ ■ ■ A (dx'" - u'" [x] dr) 

= ei-^...i^dx'^ A ■ ■ • A dx'" - nei-^^.i^u'^ [x] dr A dx'^ ■ ■ ■ A dx'" 

= n\ {dx^ A • • ■ A c/x" - c/r A d/i A • ■ ■ A . (31) 

The transition from the first to the second line here is trivial, while the detailed steps leading from 
the second to the third lines are given below, in PHjl . Even without detailed index manipulations, 
however, the form of the third line is transparent given (j^^ and the form of the second line, and 
vice versa. 

For instance, to go from the second to the third lines, first note that Q is first-order in u. Then, 
since is multi-linear and totally antisymmetric in the first partial x^^'' derivatives of the Is, the 
terms linear in must sum to give dli A ■ ■ ■ A din-i- The only thing that is perhaps not obvious 
is the relative coefficient between the two exact n- forms that appear in the third line of (jHTj) . This 
relative coefficient is carefully and directly determined below, in ()48|) . An alternate albeit indirect 
verification of the relative coefficient, as written in (|3ip. was first given in ^3] through Tahktajan's 
demonstration that the obvious extrema of the integrals of the first and second lines of (j3ip are also 
extrema of integrals of the third line. This test would fail for other relative coefficients. 

Further consideration of such integrals is warranted. We want Q to be exact^ to trivially reduce 
Jj^^j Q to Jgj^j A, i.e. to a "boundary action," that describes evolution of the d-brane. That action 
will be manifestly extremal on the solution set defined by (j27p . given the multilinear form of Q in 
terms of those first order equations of motion. The action that governs the evolution of the d-brane 
(for n = 3, see [3], and for general n, see jl4j, especially Remark 9) is 

A= n = n\ (x^ A rfx^ A • • ■ A dx" + hdr A d/2 A ■ ■ ■ A . (32) 

J Mn JdM„ 

The latter form of the action is an integral over the n — 1 dimensional world- volume swept out by the 
evolving d-brane, confirming that at any instant the d-brane is of dimension n — 2. Parameterizing 
the world-volume by coordinates ai, ■ ■ ■ , q;„-2, as well already chosen, the action is 



A = drdai ■ ■ ■ dan-2 

'dM„ 



(33) 



{n - 1)! ei^...i„x'^drx'^da^x''-' ■ --da^.^x'" 

NB In this last expression, all is are summed from 1 to n, but all js are summed from 1 to n — 1. 

The equations of motion that follow from extremizing this action are then projected versions of 
the original (jTTj) . namely 

jdrx'' - u'') d^,x'' ■ ■ ■ 9„„_,x*" = . (34) 

necessary condition for to be exact is that it be closed. This imphes the flow should be divergenceless, 
V-M= 0, as eij...i„du*i [x] A dx*^ • • • A dx'" = ei^...i^dju'^dx^ A dx*^ • . • A dx'" = i {djU^) Ei^.-.i^dx"-^ A (ix*^ • • • A dx'" . 
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This projected form of the equations of motion (|27|) allows for arbitrary reparameterizations of the 
{n — 2)-brane surface, at any given r, including r-dependent choices for the aj. Through this 
reparameterization freedom, as emphasized by Regge et al. |7j, and subsequently by Tahktajan 
[131 (especially Remark 8), there exists a r-dependent parameterization such that (jTTj) is fully 
recovered. However, in more geometrical/physical terms, the evolution of a point according to 
(j27p is fully recovered just by considering the (j34j) -driven evolution of n — 1 intersecting branes, 
appropriately configured to intersect at the point in question. This is most easily visualized for the 
string case, with n = 3. A point defined by the intersection of two appropriately chosen strings, 
with each string moving according to will itself evolve according to (jlTj). 

What is most interesting in all this is the incorporation of interactions through the use of Nambu 
brackets, particularly for the case of even dimensional phase-space d-branes for superintegrable 
systems. This elegant geometrical characterization of integrable systems with a maximal number 
of invariants has not been fully appreciated previously, I believe. Also, this novel geometrical view 
holds considerable promise for understanding the quantization of such maximally superintegrable 
systems - promise made all the more possible by recent work on the quantization of the Nambu 
bracket using both traditional Hilbert space operator and non-Abelian deformation methods of 
quantization j2]. 

Specification of the dynamics through the n-form Q does not by itself endow the brane with 
any inherent stability or cohesion. That is, the brane is tensionless so far as the above action 
dictates. The data comprising the brane evolves smoothly and remains contiguous without tearing 
only through the courtesy of the driving velocity field. If the latter is smooth enough, then so 
is the evolution of the extended data set. It may move forward along with the ambient flow and 
give the appearance of stability, much as a thin sheet of smoke, dust, or snow blown by the wind 
may stay together if the driving flow of air allows it. Moreover, for a r-independent but otherwise 
arbitrarily chosen parameterization, the d-brane responds similarly to an untethered, tensionless, 
"teflon-coated" sail (for the n = 4 case, with more challenging mental pictures for generalizations 
to higher n). Only the flow component perpendicular to the surface has any apparent effect on 
it. This would suggest that any shear in the ambient flow will cause all impacted sections of 
the surface to rotate and become parallel to u (i.e. tangential). Once parallel to u, the surface 
sections will cease being moved extrinsically by the velocity field, although intrinsic r-dependent 
reparameterizations can be used to maintain the tangential motion inherent in (j27|) . At least, that 
is the case for motion due to the action of the n-form alone. 

For classification purposes, it is sufficient to consider just the free kinetic portion of the n-form 
(Hopf term). 

uj = Ei^.-.i^dx'^ A ■ ■ • A dx*" = n\ [dx^ A • ■ ■ A rfa;") . (35) 
This alone is exact, so we obtain only a boundary action from its integral. 

/ 00 = ei^...i^x''^dx^^A- ■ ■Adx''" = (n - 1)! / drdai ■ ■ ■ dan-2Si^...i„x''^drX^^daiX''^ ■ ■ ■ da^_^x'''^ , 

(36) 

where in the last step we have again used the parameterization of the (n — l)-dimensional world- 
volume swept out by the evolving (n — 2)-brane. We have singled out one of the parameters as the 
Nambu time, r. So the free term may be used to classify the various extended structures, as for 



7 



example: 
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JM2 


Jdhh 
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JM3 


JdMs 
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f 




JdAh 


/ - 


{2N)\ 


J M2N 






{2N)\ 



J dr EijX^drX^ point particle, (37) 



c/x^ A rfp^ Adx^ Adp^ A--- A dx^ A dp^ = 2^ {2N - 1)!! / (cjs)^ 

A/2 AT A-/2]V 

/ x'AVA...Adx-Adp" maximal phase-space d-brane, 

J Af2]v-i=i9Af2]v 

where in the next to last line, UJ2 = Ylj d-^"^ ^ d,P^ is the canonical two-form on the phase-space. 

Take the simplest case of a particle in phase-space undergoing Hamiltonian flow, with n = 2. 
The phase-space velocity field in this case is given by 

u^ = ^H, u^ = -^H, ^u' + ^uP = 0. (38) 
op ox ox op 

The additional term is then 

[x] dr A dx'^ ■■■A dx^" dr A {u^dp - u^dx) = dr A {dpHdp + d^Hdx) = drAdH , (39) 
the expected exact 2-form, well-known from Hamiltonian dynamics. So the action is the usual 



A= n= xdp-pdx + 2HdT = 2 Hdr - pdx . (40) 

J M2 Jdhh Jdhh 

The method always works for the point-particle case when the flow is Hamiltonian. 
Take the second simplest case of a "vortex string" (see Regge et al. [7]), with n = 3. 

Ei^.-i^u''^ [x] dr A dx'^ ■ ■■ A dx'" ^ 2 rfr A {u'^dy A dz + vPdz A dx + u'^dx A dy) . (41) 

When the velocity field is given by a Nambu bracket, with invariant generating entries that we 
shall designate H and L, we obtain: 

ti^ = {x, H, L} = dyHdzL - d^HdyL , 
uy = {y, H, L} = d,Hd,L - d,Hd,L , 

u' = {z, H, L} = d^HdyL - dyHd^L . (42) 

For the velocity field around a vortex, these two invariants may in fact be written in terms of the 
Clebsch potentials jTUHH]. In any case, 

u^dy Adz + u^dz A dx + u^dx A dy 
= [dyHdzL — dzHdyL) dy Adz + {dzHd^L — d^HdzL) dz A dx + [d^HdyL — dyHd^L) dx A dy 
= [dxHdx + dyHdy + dzHdz) A {d^Ldx + dyLdy + dzLdz) 

= dH AdL. (43) 
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Thus the interaction with the velocity field becomes 



ei^...iji'-^ \x\ dr A dx'^ • • ■ A dx'" dr AdH AdL , (44) 

/ = 6 / dx A dy A dz — dr A dH A dL = 6 / xdy A dz — HdL A dr . (45) 

J M3 JM-i JdM'i 

Next, consider the details of the general case. We will use the generalized Kronecker symbol 
and its trace (both is and js are summed from 1 to n here). 



^n-..e,,-,. = = X 6^-:::i - 6^ x ^--.X + - (n terms), 

= ^^n-X = X C:.t • (46) 
Note from the examples the point is always to make an exact form out of 

ei,...iy' [x] dx'^ ■■■A dx'" . (47) 

But when the velocity field is given by an n-bracket, u*^ = {a;*^, /i, ■ ■ ■ , /n-i} this becomes 

ei,...i„ {x'\Ii, ■ ■ ■ , In-i} dx'^ A ■ ■ ■ A dx'" = ei,...i^ej,...j„ dj.x'^djji ■ ■ ■ djjn-idx'^ A ■ ■ ■ A dx'" 

= ei,...i„ej,...jj]l djji ■ ■ ■ djjn-idx'^ A ■ ■ ■ A dx'" 
= 6ii::i6l d,,h ■ ■ ■ djjn-idx'^ A ■ ■ ■ A dx'" 
= ^hX ^nh ■ ■ ■ d,Jn-idx'' A ■ ■ ■ A dx'" 
= {n- 1)! dijidx'^ A ■ ■ ■ A dijn-idx'" 
= (n-iy.dhA--- Adin-i . (48) 

That's all there is to it. In summary, there are two steps needed: 

ei^...i^u'^ [x] dx'^ ■ ■■ A dx'" = Ei^.-i^ /i, ■ ■ ■ , /„-i} dx''^ A • • ■ A dx'" = {n - l)\dli A ■ ■ • A . 

(49) 

Thus, for such velocity fields given by Nambu brackets, we recover (jHTj) . 

A more complete recount of the history of the ambient "fiowing fiuid" picture of the classical 
dynamics of extended objects would take note of the following. Following the seminal paper of 
Nambu, the classical theory in the general case was worked out in various stages, by Estabrook, by 
Lund, Rasetti, and Regge, by Tahktajan, and by Matsuo and Shibusa. While extended objects 
were not explicitly discussed, this picture was nevertheless initiated, formally, as a generalization 
of the usual Hamiltonian 2-form framework to higher- dimensional forms by Estabrook j3] in the 
context of Nambu mechanics 0. A physical interpretation as an extended object moving in a 
three-dimensional fiuid was given by Regge et al. [7j, for the case of vortex strings, but without 
using N-brackets. The whole business was mathematically codified for general n by Tahktajan fPS, . 
using N-brackets, but without any additional physical interpretation, and with some unnecessary 
complication involving an extended phase-space picture. More recently, Matsuo and Shibusa jH] 
have developed further the ideas of Regge et al. to higher n, and Pioline |T2] has given a summary 
of past developments, with emphasis on strings and membranes, while stressing the case of 3-forms. 
Even more recently, Curtright and Zachos |2j have discussed at length the Nambu formalism in 
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conventional phase-space contexts for superintegrable systems, have emphasized the differences 
between even and odd N-brackets, have pointed out the particularly nice properties of implementing 
Nambu mechanics with 4-brackets for models with su (2) invariance, and have carefully discussed 
the quantization of the formalism using conventional Hilbert space and/or deformation quantization 
methods. Finally, in a conference talk reviewing [2], Zachos has discussed 2-branes using 4-brackets 

m 

Only classical considerations were made above. However, the resulting mental picture evoked 
by the formalism is somewhere intermediate between the classical Hamiltonian evolution of ideal 
points, and the quantum evolution of globally well-behaved distributions (Wigner functions, say, 
defined throughout the full phase space) with support of codimension zero. Here we are led to 
envision distributions of data on surfaces of non-zero codimension, evolved quite elegantly through 
the use of a single Nambu bracket in such a way that the implicit evolution of ideal points comprising 
the surface is not necessarily just uniform time evolution, but rather evolution involving a range of 
time scales dictated by the breadth of the extended distribution in the various invariants. 

That is to say, the complete extended object evolves smoothly and uniformly under the passage 
of a single Nambu time, r. But each individual point comprising the object may evolve according 
to its own time scale such that a spectrum of conventional individual particle times t pass during 
the evolution. Aside from being at least implicit in the references cited earlier, the only similar 
classical precursor of such collectively evolved data, as far as we know, arises in the formalism 
of the phase-space surface-of-section used in the exploration of chaos and to test for integrability. 
In that formalism, the points making-up the initial surface are evolved individually, for example 
by Hamiltonian methods (see, for example, [U]). When following individual trajectories of points 
initially on the surface of section, the times required to return to the initial surface may vary between 
the various trajectories. Under Nambu time evolution, the analogue of the surface of section would 
evolve as a coherent whole for each Nambu time increment. 

The quantization of the d-brane dynamics described here is under investigation. Early results 
on quantization issues for 1-brane vortices are described in [7j. A follow-up infusion of general, more 
rigorous results for quantized 2-branes of vorticity is given in p]. Still more recent results for 1- 
branes are in |31 ^] . None of these other studies of the quantization of vorticity and 1-branes makes 
particular use of quantum Nambu brackets, however. Such usage seems far less problematic today 
than it did, say, one year ago due to the developments in When N-brackets are quantized, they 
have many surprising features, not the least of which is that the time-scale becomes dynamical in all 
but one known special case, when applied to the evolution of point-particle Wigner functions. This 
is one feature that naturally associates quantum N-brackets to the evolution of d-branes described 
above, and strongly suggests that the latter can provide a convenient bridge between classical and 
quantum behavior. A full report on this subject will have to be given elsewhere. 
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